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Восстановление разрьвньжх функций 
разрьвньюми аппроксимационнькюми сплайнами 
с использованием трапецевидньх злементов 


У статті побудований розривний апроксимаційний сплайн, коефіцієнти якого знаходяться методом 
найменших квадратів, для наближення розривної функції двох змінних з областю визначення, що 
можна розбити на прямокутні трапеції. Запропонований метод наближення можна буде використати 
для математичного моделювання розривних процесів в медичних, геологічних, космічних та інших 
дослідженнях. 
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В статье построен разрьгювньшй аппроксимационньтгй сплайн, козффициентьт которого находятся методом 
найменьших квадратов, для приближения функции двух переменньх с областью определения, которую 
можно разбить на прямоугольньтєе трапеции. Предложенньй метод можно использовать для математического 
моделирования разрьтвньїх процессов в медицинских, геологических, космических и других исследованиях. 
Ключевьге слова: разрьвная функция, разрьвньвй аппроксимационньій сплайн, 

метод найменьших квадратов, томография 


Вступ 


Задача наближення неперервних функцій неперервними сплайнами з достатньою 
повнотою описана в багатьох роботах вітчизняних та зарубіжних дослідників |1-31. 
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У вказаних роботах досліджувалося наближення неперервних функцій за допомогою 
неперервних та розривних сплайнів. Але загальної теорії таких наближень не існує. 

Задачі наближення розривних функцій виникають частіше, ніж задачі наближення 
неперервних функцій. Наприклад, в методах комп'ютерної томографії на даний час 
ніде не використовується інформація про внутрішню структуру тіла людини (шлунок 
має одну форму 1 відповідну щільність його тканин, печінка має іншу форму та іншу 
щільність тканин тощо); при дослідженні кори Землі за допомогою даних з кернів 
свердловинного буріння виникає задача відновлення внутрішньої структури кори між 
свердловинами. При цьому очевидним є той факт, що щільність грунту в різних точках 
кори є неоднорідною і найчастіше має розриви першого роду в точках поверхонь, які 
відділяють одну складову кори від іншої (чорнозем, пісок, глина, граніт тощо). Тобто 
актуальною є розробка та дослідження теорії наближення розривних функцій за допо- 
могою розривних функцій. 

У статті |4| авторами були побудовані розривні лінійні інтерполяційні сплайни для 
наближення функцій однієї змінної, що може мати розриви першого роду. У роботі |3| 
був запропонований метод наближення розривних функцій двох змінних розривними 
інтерполяційними білінійними сплайнами, а в роботі (6) - інтерлінаційними розривними 
сплайнами. Були також побудовані розривні інтерлінаційні сплайни для наближення 
функцій двох змінних, область визначення яких розбивається на прямокутні трикутни- 
ки |7|. У статті |8| запропонований метод відновлення розривних функцій двох змінних, 
з використанням трапецієподібних елементів. 

У даній роботі будуються та досліджуються апроксимаційні розривні сплайни 
для відновлення розривних функцій з областю визначення, що розбивається на прямо- 
кутні трапеції. 

Постановка задачі. Нехай задана розривна функція двох змінних /(х, у) в області 


реї0,п. Будемо вважати, що її область розбита на прямокутні. Трапеції не 
вкладаються одна в одну, а їх сторони не перетинаються. Функція /(х, у) має розриви 


першого роду на границях між цими трапеціями (не обов'язково між всіма). 
Потрібно побудувати апроксимаційний розривний сплайн для заданої розривної 
функції /(х, у). 


Побудови розривного апроксимаційного сплайна 


Якщо (х;, у;) - вузол, в якому знаходиться прямий кут прямокутника, то може 


зустрітися чотири типи трапецій (рис. 1): 
1 1 5 
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Хі-1 Хі Хіч 
хі-1 хі ЖІН 


Рисунок І - Зображення можливих трапецієподібних елементів 
з прямим кутом у вузлі (Х,, У) 


Вважаємо, що на кожній зі сторін заданих трапецій функція /(х, у) може мати 
(а може і не мати) розриви першого роду. 
Визначення. Будемо називати розривним апроксимаційним поліноміальним сплай- 


ном в області ТРІО ср,ко18 таку функцію: 
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І 
(ха 70,8 р (ха) 70, 


Раса 
І! 
ма 


(70,8 | (0 )-0), К-2 
С 20,8 б (ау), кКЕЗ 
(ха 08 дО), 4 
СОДУ -0,Уці70), Ко 
Оу а) нО уза 20), Каб 
А ураужОуїіно к- 
0) -0, Уа О кеВ 
а коефіцієнти СЯ ),ї, ) - 1,4 сплайна знаходяться методом найменших квадратів з умови 


||| СГО5 у) - 503, у) Ах з тіп 2) 
ВР, б 


Теорема. Для оператора наближення розривної функції /(х, у)є С Кн ТВО) роз- 
ривним апроксимаційним сплайном 5(х, у) вигляду (1), побудованого за допомогою 
методу найменших квадратів, на кожному елементі розбиття ТВ», і - І, т, і -1,п, спра- 


ведлива така оцінка: 
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Доведення. Проведемо доведення на прикладі трапеції ТР Тоді формула (1) 


перетвориться у такий вираз: 


ек (1) 
5(х у) У "денних а ХУ Хі У 5 уч (9 Ж сів ХО Хі і Де 5 уч (х) 
: І І 
ХОЯІ У о ре (у) ХП" Уз з 2 (У) 
Оу ОО у на ок 93 хох 
3 4) ерерчніь. ЕТО во 
5 уч (х)- р іч 5 уч (ха) - Уч о 
Розв'яжемо мінімізаційну задачу: 
хи 810 (2 
Р.(С) пе Ї 1 (стат Х- Хі у- ноз. са) Хо-Х; у- 814109 
і) І 1 
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2 
1 тай зн б зай | но РО у) | ахау -з тіп 
І І | 
з - УР Ярі злі) -Уї Хин 7 У; С 


Вираз у дужках позначимо через ЩХ, у, С). 
Випишемо систему лінійних алгебраїчних рівнянь 


др (Су дв (С дВСу дв С) 
спас а ас 


відносно невідомих С ай ке1234 


(0 
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хо У 8 1405) У; ХО Ха 
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1 1 
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1 І 1 І 
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ізамінимо 7 (х, у) інтерполяційним сплайном, побудованим на трапецієподібному еле- 
: а)» М У . 

менті ТРі; із залишковим членом К(х, у), який був виведений у роботі |3|. В резуль- 

таті отримаємо такі вирази для інтегральних членів отриманої системи, беручи до 


уваги, що зоре ахор. 


канОЮ 2 1 - 
чн б резіВ у-вун (9 (дн о КАБ Чу, Зах ах) 


йкдусннн о 


зе Хі Хі ІЗ (х) 36 
а) 2 
МИР НТ, 
хід Вч зонди у- о х) з 0 (ваза (2-2; каху нах). 
2 1 узнотж» его З 
а) 
зон ХУ Ху | (-УД0- 80109) ахау - (а ЗДА -4у) Зхачахуц) 
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2 2 і 
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а) 2 
ад бін ФО 2 1 
хна Від хох; їЕ у- 809 за 46 -4У з Б Зах; БаХуц хро Х 
хау з ; 
1 2 
хо У; Хіт 7 Хі ія ()) 12 З 


(хар ХА -- ду; ках; чЗахці). 


(2 
сову 2 1 
А хох | У ВНС ду 
177 


2 
мії І 36 
Хі Уї Хі (з19- у) 
(7 2 
М НО, 2 
о Ре У- У; ахау М п А В 
за ОО 1 9 
хо У Хі Хі 8009 г 36 
2. 


ху ВУНОО М " 
(ха 2-2 з нах; нах) 


(хо хх о хи) й і РР НЯ у 


2 (0 
хо У) (а .9 8109 - У; "8 
(1 2 
; 8 і (9 2 
Хід 8/1 хох | у- У; МОР (хар ох ДАБ- «Ку, зах; ноахца 
р а о 1. 
хо у; АЧНГЯ вл) .) 36 
Система (4) тоді буде складатися з таких рівнянь: 
4Р -4у; Зах; Захід хро Х є. оЖнтЯі 2Ь- у; ках; Бах пути 
12 3, М 6 і 
(ура АВ - Фу; За ах) (ха 0-2 з вах; КАХ і) Й 
- 72 Єна 72 Єна, 
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а) 
ал 8 09 1 
НБН вої У 8 (о 


-Г Р ко в о 
х Уї у Хі 5 р і 8 109 
хиртх; 20-2Ууз ках Бах д 4Б -4у; 4 Зах; ах хо 
З о 
Р В п за "АБ -4у; ах, заду. я 
72 і,)-1 36 і-Н1, /-1 
здо о) 
у- 81109 
Е Ї косу) з, уа; 
хі У; Хі 7 Хі 8109 
(хі -дАР-фу; зЗхачахці) (хі ар М зах, нах і) 
винос ббееЕВВнь. 
(ха "ХУ -2у з нах; Бах) 


(ху 7246 -4у з Зах; чахул) М 
езнрокненавинн цу У Ззанвио 13. бле роз оон 100 ба 


36 і 36 
ху ЗМО) дай усу» 
- ДОР посуд 
ія 841097 У; 


Хі Уї 
д (ха 0-2 з нах; ах ) (а "АБ - 4У; вах; Зах) 
ікнаньно зу Я аванс за 
(ха ХОБ-2у з нах; ках) (ху 7246 - Фу; ках; Зах) 
зн занос т 

Хі 809 хох у-у; 

- Ї | КО, у) - Ф 7 фах. 
Ще ху 81109) 


Для аналізу правих частин отриманої системи скористаємося формулою з роботи (8: 
705 у) - 505 У «0, 


де О визначається формулою (3). 
Використовуючи позначення ЦЕ | - - тах Є, збі, і» віч Єна, ун) ТА спростивши 


отримані вирази, систему (4) перепишемо у вигляді 


Са о хі) (Зхна «зу» ка дах; че ахіі) Ів | «0: (хі о хі (Зх зу» А дах; Я ах) . 
12 12 

(хар Ху ід. зах; -Захуцц) і| «0: (хар ОХ еЗуї зах; Я Захуц) 
12 12 

(Хі 7 Хо) (Зхна -Зу; ках я ахіі) 5 (ха п Ху о ЗУ з даху В аху), 
- | «0: і, 

(Са 7 Ха ЗУ Р ху даху) «о: Са о ХХ З ЗУ ах; дал) 


12 М 12 
Тобто є | «0. 
Теорема доведена. 
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Наслідок. Якщо наближувана функція /(х, у) є кусково-лінійною або кусково- 
сталою функцією в кожному трапецієподібному елементі розбиття з точками розриву 


(о, У і» і- І, т, і - І,п у випадку наближення її кусково-лінійним сплайном 5(Х, у), 


визначеним формулами (1) з невідомими С )т-1,4, Кк «1,8, що знаходяться з 
умови (2), то отримаємо точно наближувану функцію, тобто 5(Х, у) - /(Х, у), де 
ЛО», у) - А(соп5т) або 7 (х, у) з Ду Я Аух к Аз у ж Азху. 
Зауваження. Якщо 
с осі ос сід - 


2505, У, зт, у зЇп, 
або 
ебаа пн воно 


то побудований розривний апроксимаційний сплайн вигляду (1) є неперервним ліній- 
ним апроксимаційним сплайном. 
Приклад. Нехай задана функція на одиничному квадраті |0, Щ) х Ї0, 1) (рис. 2): 


хзу, 0 У«х«ї1,05«у«-0Ахчі1 
ода 15-4х-у?,  0«х« 0,5, 05«у«04хн 06. 

0.5, 0 « х« 0,5, 04-04Ах « у « 0,5 

-х31.5, 0 5«х«1,0Ах«у« 0.5 


Тобто на лініях фігури, зображеної на рис. 2 а), функція /(х, у) має розриви пер- 
шого роду. Нехай задані лінії: 
х-0,х, 20.5, х3 -1, 
2 23019 2 2 2 
У з, оон ре. го Я 
Вони розбивають область визначення функції /(х, у) на вісім трапецієподібних 
елементів (один з кутів обов'язково є прямокутним). 


с 


МИР 


а) 6) 
Рисунок 2 - Графічне зображення: а) області визначення функції /(х, у); 
6) функції /(х, у) 
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Спочатку побудуємо розривний інтерполяційний сплайн на заданій сітці, який, 


наприклад, для трапеції ТР) задається формулою 


з РО пу «фев 
505 уз РОор ун вні, ЧО 77 пн а: у-вунод - 
Пет (1) іні» 7) ря (0 
у-у; хохо | УЗУ; | 
М РАЗ А РУ ЕН о о он внбіні ес 
/ (1) / 
ХХ 84 09- У; -х вод У 


Отримаємо інтерполяційний сплайн, графічний вигляд якого наведений на рис. 3. 


Рисунок 3 - Графічний вигляд розривного інтерполяційного 
сплайна (світлий колір) та заданої функції (темний колір) 


Далі визначимо максимальне відхилення наближуваної функції /(х, у) від побу- 
дованого сплайна 5(Х, у): 


тах| (5, У - 5( У)| я 0.3. 


Тепер побудуємо розривний апроксимаційний сплайн на заданій сітці, який, 


наприклад, для трапеції ТВО задається формулою 


5 (х, уз фрі-чн о би ; чн у-8 ще б Х- Х; у-8 но 
1 
Ро кана У У; -Со-яая У У 


4 1 
ходнівлОд-у || чат в9д(-У; 


Після знаходження коефіцієнтів цього сплайна з умови (2), отримаємо сплайн, 
наведений на рис. 4. 
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Рисунок 4 - Графічний вигляд розривного апроксимаційного 
сплайна (світлий колір) та заданої функції (темний колір) 


Як на рис. 3, так і на рис. 4 бачимо, що задана функція наближується точно там, 
де вона задана константою або лінійною функцією (графіки функції /(х, у) та набли- 


жуючих сплайнів точно збігаються). 
Далі визначимо максимальне відхилення наближуваної функції /(х, у) від побу- 


дованого сплайна 5(х, у): 
тах |/ (х, у) - 5 (х, у)| я. 0.08. 


Як бачимо, побудований розривний апроксимаційний сплайн наближує розривну 
функцію краще, ніж інтерполяційний. Побудовані розривні сплайни точно наближують 
ту частину функції, де вона є постійною або лінійною, що 1 підтверджує викладену 
вище теорію. 


Висновки 


У роботі пропонується метод побудови розривного апроксимаційного сплайна для 
наближення функції з розривами першого роду та область визначення яких розбита на 
прямокутні трапеції. Причому побудовані розривні сплайни включають в себе як частин- 
ний випадок класичні неперервні сплайни першого степеня на заданій сітці вузлів. 

В подальшому авторами планується знайти оцінку похибки наближення побудо- 
ваної розривної конструкції та застосувати розроблені методи до відновлення фантома 
Шепа-Логана в комп'ютерній томографії. 
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О.М. Іуїууп, У.І. Рег5Ніпа 
Кезіотатоп ої Дізсопіїписоиз Кипспоп5 Бу Різсопіїписиз 


Арргохітаїіопа! Уріїпез ул С5е ої Ткарегез Еіетепіз 

І Фе уогК, Фе ргобіет ої арргоасі ої даїзсопійпиоця Гапсоп ої бо уапіабіе8 іп Фе агеа 
гергезепійпя а шпії 5дцаге 15 сопз5ідегей. Її 15 5заррозед, Шаг Ше агеа БгеаК5 їпіо гесіапеціаг 
ігарегез8, місії аге пої риї еасп оШег, апа еїг рагіїез аге пої сго55ед. БВипсйоп Па8 гаріиге5 ої 
ре Пг8с 5огі оп Богаег5 Бебууееп ШПезе ігарегез (пої песез5агіїу Бебмееп а). п ре агісіє, її раз 
Бееп соп5ігисіва Фізсопійпиоц5 арргохітанопаї! роїупопаа! 5ріїпе, мурісб соеййсіепія аге Гоцпа 
Бу Фе теоа ої Фе Іеахі 5дпагез, Їог арргоасії ої Фе зеї Ффізсопійпибсицз Гипсбоп. Тпе Шеогет 
ої ап езйтайоп ої Фе сопяїгисітеа фізсопіписи5 дезієп оп еасп еЇетлепі ої зрійштє 15 ргоуей. 
І Фе могк, сопайоп5, мПпісп арргохітацопа! дїзсопйпиоц5 5ріїпе 15 сопіїпиоця аї, аге 
тезийеа, Фа 15 Фе яепегаї їреогу ої арргоасП ої ре 5еї Гипсйоп5 у/ПпісП 5ресіа! сазе 15 а! Фе 
Кпомп ШФеогу ої арргоасп Бу сіазкіса| сопійпиоц5 зріїпе5 15 деуеїоред. Те ехатріе ої 
арргоасі Фе Типсбоп ої ємо уагіабіе8 Бу іпегроіацопа! а43зсопійпиоия 5ріїпе раз Бееп деуеіоред 
Бу айфог5 Беїоге, апа сопяігиства їп ре агісіе арргохітайопаї! д4ізсопипиосиця зрітпе 18 гезиітей. 
Тре гесеїуеа гезиія ої арргоасі аге аїзо апаїугей. ТРре 5ег ехатріе аї5о сопіїття Фе ШФеогу 
згакед їп агіісіє. 

Тре деуеіоред тефоа ої арргоасії аззитез, Фас Типспоп гиріагез аге кпом/п, апа Шеу 
соїпсіде ул Кпої5 аїзсопійписицз арргохітайопаї зріїпе. Те ашфогз ріап Іакег іо демеїор 
теШодзя ої арргоасп ої азсопійпиосиця Гипсйопя її шпкпомуп гиріиге8з, апа іо арріу Ше 
деуеїоред теодя іо Ше Чесізіоп ої а (муо-дітепзіопаї! ргобіеті сопіриїег ої гоптоєтарбу, 
патеїу, їо гез'огайоп ої Кпомуп 5рерр-І.оєап ррапіот. ТПе оНегей арргоасі сап Бе цзеа Їог 
тафетайсаї подеШтпе, ої афФізсопіпиоця ргосе85е5 їп пледїса!, єеоіовісаї, зрасе апа офег 
гезеагсрез. 
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